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The Dynamical Cluster Approximation (DCA) is used 
to study non-local corrections to the dynamical mean field 
phase diagram of the two-dimensional Hubbard model. 
Regions of antiferromagnetic, d-wave superconducting, 
pseudo-gapped non-Fermi liquid, and Fermi liquid behav- 
iors are found, in rough agreement with the generic phase 
diagram of the cuprates. The non-local fluctuations beyond 
the mean field both suppress the antiferromagnetism and 
mediate the superconductivity. 

Introduction The rich phenomenology of high-T c 
superconductors has stimulated strong experimen- 
tal and theoretical interest in the field of strongly cor- 
related electron systems. Common to all high-T c sys- 
tems is the presence of antiferromagnetic ordering in 
undoped samples in proximity to a superconducting 
phase with a d-wave order parameter and the normal 
state pseudogap dominating the physics in underdoped 
samples. A successful theory must describe all these 
fundamental features at the same time. 

The 2D Hubbard model in the intermediate coupling 
regime or closely related models like the t-J model 
are believed to capture the essential physics of the 
high-r c cuprates 0. The antiferromagnetic phase of 
the cuprates is well understood. In the strong cou- 
pling limit U ^ W, where U is the Coulomb repul- 
sion and W the bare bandwidth, the undoped Hub- 
bard model reduces to the Heisenberg model, which 
has been proven to describe the low energy spin fluctu- 
ations of the cuprate parent compounds. However, off 
half-filling there is no complete understanding of the 
superconducting phase or the normal state pseudogap 
in the intermediate coupling 2D Hubbard model. 

Finite size quantum Monte Carlo (QMC) calcula- 
tions for the doped 2D Hubbard model in the interme- 
diate coupling regime with U < W, support the idea of 
a spin fluctuation driven interaction mediating d-wave 
superconductivity ||. However, the fermion sign prob- 
lem and the fact that the number of degrees of freedom 
grows rapidly with the lattice size, limits these calcu- 
lations to temperatures too high to study a possible 
transition j|. These calculations are also restricted 
to relatively small system sizes making statements for 
the thermodynamic limit problematic, and inhibiting 
studies of the low energy physics. 

These shortcomings do not apply to the Dynamical 
Mean Field Approximation (DMFA), which is by con- 
struction in the thermodynamic limit. Unfortunately, 



the lack of non-local dynamics in the DMFA inhibits a 
possible transition to a state with a non-local (d-wave) 
order parameter. 

However, it is well known within phenomenologi- 
cal theories that short-ranged antiferromagnetic spin 
fluctuations mediate pairing with d-wave symmetry 
and cause a pseudogap in underdoped samples 
There is experimental evidence that the correlation 
length of dynamical spin fluctuations in the optimally 
doped cuprates is very small Therefore, micro- 
scopic studies, which account for short-ranged dynam- 
ical correlations in addition to the local correlations of 
the DMFA are relevant and might succeed in describing 
the physics of the cuprates. 

The recently developed Dynamical Cluster Approx- 
imation (DCA) @-§ is a fully causal approach which 
systematically incorporates these non-local corrections 
to the DMFA by mapping the lattice problem onto an 
embedded periodic cluster of size N c . For N c = 1 the 
DCA is equivalent to the DMFA and by increasing N c 
the dynamic correlation length can be gradually in- 
creased while the calculation remains in the thermody- 
namic limit. Previous DCA calculations have indicated 
the presence of an extended d-wave superconducting 
phase at moderate doping. [|§|-|lof 

In this manuscript we present calculations of the full 
phase diagram of the 2D Hubbard model studied with 
the DCA. To solve the cluster problem we use QMC. 
We choose iV c = 4, the smallest cluster that includes 
non-local corrections while preserving the full transla- 
tional and point group symmetries of the lattice. The 
results are compared to DMFA calculations, N c = 1, to 
illustrate the effect of the initial non-local corrections. 

Formalism A detailed discussion of the DCA for- 
malism was given in previous publications |^-|s|JTl|. 
The compact part of the free energy is coarse-grained in 
reciprocal space, projecting it onto a finite-sized clus- 
ter of iV c points embedded in a self-consistently deter- 
mined host. The cluster problem is solved using the 
Hirsch-Fye impurity algorithm fl2]| modified to simu- 
late an embedded cluster [jll] , and the spectra are ana- 
lytically continued with the maximum entropy method 
[ [l3[ . Once the cluster problem has been solved, lattice 
susceptibilities may be calculated All calculations 
are done in the normal, paramagnetic state. We search 
for continuous phase transitions indicated by the diver- 
gence of the corresponding susceptibilities. 
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FIG. 1. Inverse antiferromagnetic susceptibility versus 
temperature for U — 2. The lines are fits to the function 
1/Xaf{T) = b(T-T N y. For N c = 1 7 « 1, the mean-field 
value. For N c = 4 non-local fluctuations suppress the tran- 
sition, so that 7 increases and Tjv decreases (see inset). 

The Hubbard model is characterized by a near- 
neighbor hopping t and an one-site repulsion U . We 
choose t — 1/4 to establish a unit of energy and choose 
U = W = 2 which is sufficiently large that for iV c > 4, 
a Mott gap is present in the half-filled model. The 
phase diagram of the Hubbard model, calculated with 
DMFA, displays a range of behaviors including Mott 
insulating, antiferromagnetic and Fermi liquid regimes 
fl4| . The inclusion of non-local corrections yields sig- 
nificant changes to the phase diagram, including the 
enhancement of the Mott phase at half filling |l(| , the 
suppression of antiferromagnetism, the introduction of 
a d-wave superconducting phase and non-Fermi liquid 
behavior. These different regimes are delineated by cal- 
culating the antiferromagnetic (Fig. [l]), bulk magnetic 
(Fig. ||) and d-wave superconducting (Fig. ||) suscepti- 
bilities, and the single-particle self energy (Fig. |J). 

Antiferromagnetism In Fig. [j], the inverse antifer- 
romagnetic susceptibility is plotted versus temperature 
for dopings 5 = and 0.05 and N c = A and 1. We plot 
the inverse susceptibility to show that in contrast to 
finite-sized system calculations, the susceptibility di- 
verges for low doping at T = TV- This indicates a tran- 
sition to an antiferromagnetic phase. In the DMFA for 
the 2D model, or as we found previously for the infinite- 
dimensional model fl4|| , the antiferromagnetism per- 
sists to relatively high temperatures and dopings. The 
non-local dynamical fluctuations, included in the DCA 
for N c > 1, strongly suppress the antiferromagnetism. 
Their effect becomes pronounced for low temperatures 
and dopings. For example, when 5 = 0, the -/V c = 1 and 
N c = 4 AF susceptibilities are identical at high temper- 
atures due to the lack of non-local correlations, but sep- 



arate as the temperature is lowered. The N c = 1 sus- 
ceptibility diverges with mean-field exponent of about 
one; whereas the N c = 4 result diverges at a much 
lower temperature with a larger exponent. Consistent 
with the Mermin Wagner theorem, Tn continues to fall 
for large values of iV c (not shown). 




FIG. 2. The single-particle density of states N(u>) for 
U = 2, T = 0.023 and N c = 4. The inset shows the bulk 
susceptibility as a function of temperature. For 5 < 0.2 a 
peak develops at T = T* accompanied by the evolution of 
a pseudogap in the DOS for T < T* . 

The Pseudogap and Non-Fermi Liquid Behavior 
The bulk (k = 0) magnetic susceptibility and single- 
particle density of states (DOS) display evidence of a 
pseudogap for low doping 5 < 0.2. We show the bulk 
magnetic susceptibility in the inset to Fig. || for three 
different dopings. For low to intermediate doping, it 
develops a peak at low temperatures, defining a tem- 
perature T*. T* < T N {N C = 1), the mean field tran- 
sition temperature (see Figs. || and |l|). For N c > 1, it 
defines the temperature where short-ranged spin cor- 
relations first emerge. The underdoped bulk suscepti- 
bility data, 5 < 0.075, may be scaled onto one curve by 
plotting versus T/T* (not shown). A similar peak or 
downturn and scaling is seen in the Knight-shift data 
of the cuprates (I) . The downturn of the susceptibility 
is accompanied by a loss of states near the Fermi en- 
ergy. For temperatures T < T* , a pseudogap begins to 
develop in the DOS, as shown in Fig. |[ The pseudogap 
is widest, measured from peak to peak, at low doping, 
and vanishes for 5 > 0.2. The depth of the pseudogap 
is greatest when 5 m 0.05, and it vanishes as 5 — > 0, 
where it is replaced by a Mott gap of width w U. 
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FIG. 3. The single-particle self energies at the DCA 
cluster momenta, plotted versus frequency for U = 2. 
T = 0.023, N c = 4 and 5 = 0.05 (left) and 5 = 0.25 (right). 
The self energy at (7r,0) changes from non-Fermi liquid like 
at doping 5 = 0.05 to Fermi liquid like at S = 0.25. 

The slow fall of T* with doping, indicates that the 
short-ranged spin correlations diminish slowly upon 
doping. Thus, even in the region T* > T > TV the 
antiferromagnetic correlations still have a significant 
effect. This is supported by the behavior of the self 
energy for temperatures T < T* . Once the DCA algo- 
rithm is converged, the lattice self energy is calculated 
by interpolating the cluster result on to the full lattice 
Brillouin zone. Thus the lattice self energy at any k 
is dominated by the cluster self energy at the nearest 
cluster momentum. For a Fermi liquid, the self en- 
ergy S(k,o;) ~ (1 — 1/Z)uj — ibuj 2 where b > and 
l/Z > 1. Our results show that near half filling, the 
self energy displays non- Fermi liquid behavior. This is 
illustrated in Fig |3J where we plot the low- frequency self 
energy at the DCA cluster momenta for S = 0.05 and 
T = 0.023. For momentum points near k = (7r,0), the 
imaginary part of the self energy crosses the Fermi en- 
ergy almost linearly. Concomitant with this behavior 
is a pseudogap of width « |J| ss 4i 2 /£7 in the single- 
particle spectra A(k, uj) for momenta near k = (7r, 0) 
(not shown). 

The pseudogap and the anomalies in the self energy 
vanish when T* falls to zero. Here, as shown on the 
right of Fig. ^, the self energy becomes Fermi-liquid- 
like with quasiparticle weight Z « 1/2. A systematic 
study of the evolution of the single-particle spectra and 
the Fermi surface will be presented elsewhere. 

It is important to stress that the pseudogap, the 
downturn of the bulk magnetic susceptibility and the 
non- Fermi liquid behavior in the self energy are absent 
when N c — 1 due to the lack of non-local fluctuations. 

Superconductivity We searched for many different 
types of superconductivity, including s, extended-s, p 
and d-wave, of both odd and even frequency and we 
looked for pairing at both the zone center and cor- 
ner. Of these, only the odd-frequency s-wave and even- 
frequency d-wave pair field susceptibilities at the zone 



center were strongly enhanced, and only the d-wave 
susceptibility diverged. This is illustrated in Fig. [| 
where the odd-frequency s-wave and even frequency 
d-wave q = susceptibilities are plotted versus tem- 
perature for U — 2 and S — 0.05. The s-wave and 
extended s-wave q = even frequency susceptibilities 
are also plotted for comparison. 




FIG. 4. The s-wave, extended s-wave, and d-wave even 
frequency and the odd-frequency s-wave q = suscepti- 
bilities versus temperature for U = 2.0, S — 0.05, and 
N c = 4. Pairing is found only in the even-frequency 
q = d-wave channel. In the inset the inverse d-wave 
pair-field susceptibility is plotted versus temperature for 
two different dopings and cluster sizes. The line is a fit to 
l/P d (T) = b(T - T c )"< with T c = 0.021 and 7 = 0.72. 

The behavior of the d-wave pair-field susceptibility 
as a function of temperature for N c — 1 and 4 and 
(5 = and 0.05 is shown in the inset to Fig. || For 
N c = 1 there is no tendency towards pairing. For the 
DMFA there is no pairing with symmetries lower than 
the lattice symmetry (i.e., p-, d-wave, etc.) p5| . 

d-wave pairing is strongly enhanced for N c = 4 over 
the corresponding DMFA results. However, for <5 = 
the inverse susceptibility rises abruptly as the temper- 
ature is lowered and the Mott gap opens in the DOS. 
The Mott gap becomes more pronounced as N c in- 
creases |l6|l , so that for larger clusters the gap pre- 
vents superconductivity even for U < W. If charge 
excitations are gapped, then pairing is suppressed. At 
half filling, for U = 2 the gap is of order U, and 
is much larger than the magnetic exchange energy 
I J I - At 2 /U = 0.125, so that the opening of the Mott 
gap will suppress any magnetically mediated pairing. 
Away from half filling the width of the pseudogap in 
the charge excitation spectrum is much smaller, on the 
order of J, so magnetically mediated pairing is possi- 
ble. For N c = 4 and 5 = 0.05, the d-wave pair field 
susceptibility diverges at T c as 0.021, with an exponent 
which is less than one, indicating that the fluctuations 
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beyond DMFA which suppress the antifcrromagnetism 
are also responsible for pairing. 
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romagnetism, the emergence of non- Fermi liquid (pseu- 
dogap) and d-wave superconducting regimes. The crit- 
ical exponent of the pair-field susceptibility is smaller 
than one, whereas the antiferromagnetic susceptibil- 
ity diverges with a critical exponent larger than one. 
This indicates that the same fluctuations that suppress 
antifcrromagnetism upon doping, mediate pairing. At 
half-filling the formation of the Mott gap of width J 
suppresses pairing. 
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FIG. 5. The temperature-doping phase diagram of the 
2D Hubbard model calculated with QMC and DCA for 
N c = 4, U = 2. Tjv and T c were calculated from the diver- 
gences of the antiferromagnetic and d-wave susceptibilities, 
respectively. T* was calculated from the peak of the bulk 
magnetic susceptibility. 

The phase diagram of the system is shown in Fig. [|. 
In addition to the d-wave superconducting and antifer- 
romagnetic phase boundaries, we include T* , the pseu- 
dogap temperature fixed by the the peak bulk suscep- 
tibility. At low temperatures, it serves as a boundary 
separating the observed Fermi liquid and non-Fermi 
liquid behavior. For T < T* and S < 0.2 the self 
energy shows non-Fermi liquid character for the parts 
of the Fermi surface closest to k = (w, 0); whereas, 
the low-temperature self energy is Fermi-liquid like for 
S > 0.2. The d-wave transition temperature is maxi- 
mum at 5 ks 0.05. The superconductivity persists to 
large doping, with T c dropping very slowly. In contrast 
to experimental findings, the pairing instability (pre- 
ceded by an AF instability) persists down to very low 
doping. One possible reason for this is that the model 
remains very compressible down to very low doping 
S ~ 0.025. This could be due to the lack of long-ranged 
dynamical spin correlations or stripe formation which 
could become more relevant as N c increases or when 
multiple Hubbard planes are coupled together. The ef- 
fect of such additional non-local corrections (N c > 4) 
is presently unknown. However, we believe that a fi- 
nite mean-field coupling between Hubbard planes will 
stabilize the character of the phase diagram presented 
here as N c increases. 

Conclusion We have used QMC and DCA to study 
the effect of the initial non-local corrections on the 
phase diagram of the 2D Hubbard model at intermedi- 
ate coupling U = W. The corrections make significant 
changes, including a strong suppression of the antifer- 
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